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Abstract 

The theory of symmetries of systems of coupled, ordinary differential equations 
(ODE) is used to develop a concise algorithm in order to obtain the entire space of 
solutions to vacuum Bianchi Einsteins Field Equations (EFEs). The symmetries 
used are the well known automorphisms of the Lie algebra for the corresponding 
isometry group of each Bianchi Type, as well as the scaling and the time re- 
parametrization symmetry. The application of the method to Type Vllh results 
in (a) obtaining the general solution of Type VIIq with the aid of the third Painleve 
transcendental Pur, (b) obtaining the general solution of Type Vllh with the aid 
of the sixth Painleve transcendental Pyi i( c ) the recovery of all known solutions 
(six in total) without a prior assumption of any extra symmetry; (d) The discovery 
of a new solution ( the line element given in closed form) with a G3 isometry group 
acting on T3, i.e. on time-like hyper-surfaces, along with the emergence of the line 
element describing the flat vacuum Type VIIq Bianchi Cosmology. 
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1 Introduction 



The idea of using the group of automorphisms in order to have a unified development of 
Bianchi Cosmologies has a long history [1] . In that direction Harvey [2] was the first who 
found the automorphisms of all three-dimensional Lie Algebras, while the corresponding 
results for the four-dimensional Lie Algebras have been reported in [3]. Jantzens tangent 
space approach sees the automorphism matrices as the means for achieving a convenient 
parametrization of a full scale factor matrix in terms of a desired, diagonal matrix [I]. 
Samuel and Ashtekar were the first to look upon automorphisms from a space viewpoint 
[5]. The notion of Time-Dependent Automorphism Inducing Diffeomorphisms (A.I.D.'s), 
i.e., coordinate transformations mixing space and time in the new spatial coordinates 
and inducing automorphic motions on the scale-factor matrix, the lapse, and the shift 
has been developed in [6] . The use of these covariances enables one to set the shift vector 
to zero without destroying manifest spatial homogeneity. At this stage one can use the 
"rigid" automorphisms, i.e. the remaining "gauge" symmetry, as Lie- Point Symmetries 
of the EFE's in order to reduce the order of these equations and ultimately completely 
integrate them [7j. The present work of ours consists in the application of this method 
to the case of vacuum Bianchi Type Vllh Cosmology. The method is recapitulated in 
section 2 while its application to the above mentioned type, resulting in the exhaustive 
discovery of the entire solution space, is given in section 3. In section 4 we discuss our 
results and give a brief description of the solution space in the form of two tables. 

2 The Method 

As it is well known, for spatially homogeneous space-times with a simply transitive action 
of the corresponding isometry group 0,0, the line element, assumes the form 

ds 2 = (N a N a — N 2 ) dt 2 + 2 N a <r a i dx l dt + 'y a p<T a i <T P j dx i dx j (2. 1) 

where the 1-forms erf, are defined from: 

do" = C^y A,^ a\ 3 - a% = 2 C%a\a P r (2.2) 

Then the field equations are (see e.g. 0): 

E = K a/3 K a p -K 2 -R = (2.3) 

E a = K» a C% e - K\C* m = (2.4) 

E a p = K aP + N (2K\K T(i - KK a/3 ) + 2N P (K av C\ + Kp v C\ p ) - iVR Q/3 = (2.5) 
where 

Kap = ~ (j a p + 2 lav C\N" + 2 1Pv C\ p NP) (2.6) 
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is the extrinsic curvature and 

R a/3 = C\ T C\ vlaKl p Xl ° v l T » + 2 C^ X C\ K + 2 C» aK C» pxl , u r X + 

(2.7) 

2 C\ K C^ vlaXl ™ + 2 C\ K C^ vlpxl ™ 

the Ricci tensor of the hyper-surface. 

In [B] particular space-time coordinate transformations have been found, which re- 
veal as symmetries of (12.31) . (12 .4ft . (12. 5ft the following transformations of the dependent 
variables N, N a , : 

N = N,N a = A P Q (N p + lpa P CT ), % v = K\ hP v laP (2.8) 

where the matrix A and the triplet P a must satisfy: 

= C\ v k%k\ (2.9) 



2P»C\ v k% = k a p (2.10) 

For all Bianchi Types, this system of equations admits solutions which contain three 
arbitrary functions of time plus several constants depending on the Automorphism group 
of each type. The three functions of time, are distributed among A and P (which also 
contains derivatives of these functions). So one can use this freedom either to simplify 
the form of the scale factor matrix or to set the shift vector to zero. The second action 
can always be taken, since, for every Bianchi type, all three functions appear in P a . 

In this work we adopt the latter point of view. When the shift has been set to zero, 
there is still a remaining "gauge" freedom consisting of all constant A^ (Automorphism 
group matrices). Indeed the system (12.91) . (12.101) accepts the solution A^ = constant, 
P a = 0. The generators of the corresponding motions = A a p A^„ 7^/3 , induced in the 
space of dependent variables spanned by Jap's (the lapse is given in terms of 7^/3, 
by algebraically solving the quadratic constraint equation ), are [IB] : 

x m = y ma %,^ (2.11) 

with A satisfying: 

Now, these generators define a Lie algebra and each one of them induces, through 
its integral curves, a transformation on the configuration space spanned by the 7 aj a's. If 
a generator is brought to its normal form (e.g. ^-), then the Einstein equations, writ- 
ten in terms of the new dependent variables, will not explicitly involve 2j. They thus 
become a first order system in the function Zi [TT] . If the above Lie algebra happens 
to be abelian, then all generators can be brought, to their normal form simultaneously. 
If this is not the case, we can diagonalize in one step the generators corresponding to 
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any eventual abelian subgroup. The rest of the generators (not brought in their normal 
form) continue to define a symmetry of the reduced system of EFE's if the algebra of 
the -Xm's is solvable [TJ]. One can thus repeat the previous step, by choosing one of 
these remaining generators. This choice will of course depend upon the simplifications 
brought to the system at the previous level. Finally if the algebra does not contain any 
abelian subgroup, one can always choose one of the generators, bring it to its normal 
form, reduce the system and search for its symmetries (if there are any). Lastly, two 
further symmetries of (12.31) . (12.41) . (12.51) are also present and can be used in conjunction 
with the constant automorphisms: The time reparameterization t —>■ f(t) + a, owing to 
the non-explicit appearance of time in these equations, and the scaling by a constant 
la/3 ^lap as can be straightforwardly verified. Their corresponding generators are: 

(2.13) 
(2.14) 

well as with the -Xvn's, as it can be 
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These generators commute among themselves, as 
easily checked. 



3 Application to Bianchi Type VI Ih 

We are now going to apply the Method, previously discussed, to the case of Bianchi 
Type Vllh- For this type the structures constants are 

° 13 ~ — ° 31 — ° 23 — — ° 32 — ~ n 

C 32 = ~~ C 23 = & 13 = ~~ & 31 = 1 (3-1) 

C a p = for all other values of a[3^ 

Using these values in the defining relation (12. 2p of the 1-forms af we obtain 

(0 e hx sinx e hx cosx \ 
e hx cosx —e hx sinx (3.2) 
\ J 

The corresponding vector fields (satisfying [£ a , ^] = \ C 1 a ^ 1 ) with respect to which 
the Lie Derivative of the above 1-forms is zero are: 

£i = d y & = d z £ 3 = d x + (z - h y)d y -(y + h z)d z (3.3) 

The Time Depended A.I.D.'s are described by 

/ ce hP ^ cos P(t) ce hP ®sinP(t) x(tf 
A a p = i-ce hP ®mnP(t) ce hP ^cosP(t) y(t) ] (3.4) 
\ 1 
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and 



x(t)P(t) + h 2 x(t)P(t) - h x(t) + y(t) 
2 (1 + h 2 ) ' 

y(t)P(t) + h 2 y(t)P(t)-hy(t)-i(t) P(t) 
2(l + h 2 ) '2 



(3.5) 
(3.6) 



where P(t),x(t) and y(t) are arbitrary functions of time. As we have already remarked 
the three arbitrary functions appear in P a and thus can be used to set the shift vector 
to zero. 

The remaining symmetry of the EFE's is, consequently, described by the constant 
matrix: 

/ e Sl s 2 s 3 \ 
M = 1-S2 e Sl s 4 (3.7) 
V \) 

where the parametrization has been chosen so that the matrix becomes identity for the 
zero value of all parameters. 

Thus the induced transformation on the scale factor matrix is j a p = M^Mg^y^, 
which explicitly reads: 

' 7n = e 2si 711 - 2 e Sl s 2 712 + s\ 722 
712 = e 2si 712 - s 2 2 712 + e Sl s 2 (711 - 722) 



(3.* 



7i3 = e Sl (s 3 711 + s 4 712 + 713) - s 2 (s 3 712 + S4 722 + 723) 

722 = e 2si 722 + 2 e Sl s 2 712 + s 2 2 711 

723 = e Sl (s 3 712 + s 4 722 + 723) + s 2 (s 3 711 + s 4 712 + 713) 

k 733 = s 3 2 7xi + 2 s 3 (s 4 712 + 713) + s 4 2 722 + 2 S 4 723 + 733 

The previous equations, define a group of transformations G r of dimension r = 
dim(Aut(VIIh)) = 4. The four generators of the group, can be evaluated from the 
relation: 

8 (3.9) 

(3.10) 
(3.11) 



where A = {1,2,3,4}. Applying this definition to f)3.8p we have the generators: 

9 n d d n d d 

Xi = 2 711- h 2 712- h 7137^ h 2 722 h 723 — 

0711 0*712 0713 722 723 
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X 2 = -27127^ h (711 - 722) t; — 

0711 9712 



d d d 

723^ h 2 712^ h 7i3^ 

0713 0722 0723 



<9 <9 d 
X 3 = 7li t; V 1\2T\ V 2 7i3 



^7l 3 ^723 #733 

d d d 

ll2T{ h 722 7; h 2 723" 



(3.12) 
(3.13) 



0713 0723 0733 

The algebra g r that corresponds to the group G r has the following table of commu- 
tators: 

[X 1 , X 2 ] = 0, [Xi, X 3 ] = X 3 , [Xx, X 4 ] = X 4 , ,g 
[X 2 ,X 3 ] = — X 4 , [X 2 ,X 4 ] = X 3 , [X 3 ,X 4 ] =0 

As it is evident from the above commutators ( 13. 14ft the group is non-abelian, so we 
cannot diagonalize at the same time all the generators. However, if we calculate the 
derived algebra of g r , we have 



g r , = {[X A ,X B ] : X A ,X B E g r } => g r , = {X 3 ,X 4 } 
and furthermore, it's second derived algebra reads: 

g r „ = {[X A ,X B ] : X A ,X B Eg r ,}^ g r „ = {0} 



(3.15) 



(3.16) 



Thus, the group G r is solvable since the g r „ is zero. As it is evident X 3 ,X 4 ,Y"2 
generate an Abelian subgroup, and we can, therefore, bring them to their normal form 
simultaneously. The appropriate transformation of the dependent variables is: 



' 7ll = e Ul ~ U(i 



712 = e Ul m u 2 



7i3 = e 



m—ue 



(u 3 + U 2 U 5 ) 



(3.17) 



722 = e Ul U6 m 4 



723 = e Ul u& {u 2 u 3 + u 4 u 5 ) 



" < - U6 (e" 6 + u 3 2 + 2 u 2 u 3 u b + u A ul) 



y 733 = e 

In these coordinates the generators Y 2 ,X A assume the form 



X 2 



x 1 



dui 



x 4 



d 

<~)(ir, 



x 3 



a 

du 3 



(l + 2ul-u 4 )£ 2 -u 5 £ 3 +2(u 2 + u 2 u 4 )£ I +u 3 £ o +2u 2 £- ti (3.18) 
— u - — u - — 2 - 

* dus " 8115 du(j 

Evidently, a first look at (13.171) gives the feeling that it would be hopeless even to 
write down the Einstein equations. However, the simple form of the first three of the 
generators ( 13.181) ensures us that these equations will be of first order in the functions 
Ui, U3 and U5. 
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3.1 Description of the Solution Space 

Before we begin solving the Einstein equations, a few comments on the allowable range 
of values for the functions Ui,i = 1, . . . , 6 will prove very useful. 
The determinant of j a p, is 

detfr a p} = e 3ui ~ 2ue {-u\ + u 4 ) (3.19) 

so we must have U4 > u\ . 

The two linear constraint equations, written in the new variables (13.171) . give 

E l = => ^ e- U6 ((3 h - u 2 ) u 3 + (3hu 2 - u 4 ) Us) = (3.20) 
E 2 = 0^ i e -" 6 ((1 + 3 hu 2 ) u 3 + (u 2 + 3 hu 4 ) u 5 ) =0 (3.21) 

This system admits only the trivial solution, since the determinant of the 2x2 matrix 
formed by the coefficients of u 3l u^ becomes zero only for the forbidden value M4 = u 2 . 
We thus have 

U3 = h, u b = k 5 (3.22) 

Now, these values of ^3,^5 make 713,723 functionally dependent upon 711,712,722 (see 
(I3.17P ). It is thus possible to set these two components to zero by means of an appropriate 
constant automorphism. 

We therefore can, without loss of generality, start our investigation of the solution 
space for Type Vllh vacuum Bianchi Cosmology from a block-diagonal form of the scale- 
factor matrix (and, of course, zero shift) 

f ln 712 \ 

la/3 = I 712 722 (3.23) 

733/ 

These unknown functions of time have to satisfy the quadratic and the third linear 
constraint, as well as the spatial EFE's. As we have earlier remarked, since the algebra 
(I3.14p is solvable, the remaining (reduced) generators X±,X 2 (corresponding to block- 
diagonal constant automorphisms) as well as Y 2 continue to define a Lie-Point symmetry 
of the reduced EFE's and can thus be used for further integration of this system of 
equations. 

The remaining (reduced) automorphism generators are 

X 1 = 2 711 — + 2 7i 2 — - + 2 722 — (3.24) 

»7ll O712 C?722 

d . d d 
X 2 = -2 712 h (711 - 722) ^ h 2 T12 



d7n ^712 oj. 



22 
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The appropriate change of dependent variables which brings these generators -along with 
Y 2 - into normal form, is described by the following scale-factor matrix : 

'i e ui+2i*i (i -2u 2 sin2u 4 ) e Ul+2u6 u 2 cos2w 4 

y a p= ( e ni+2u6 u 2 cos2w 4 5 e Ul+2 " 6 (1 + 2 u 2 sin 2 u 4 ) | (3.25) 

e' 

The generators are now reduced to 

d d d 

Y 2 = - — , X 2 = —,X 1 = -— (3.26) 

OU\ OU 4 OUq 

indicating that the system will be of first order in the derivatives of these variables. 
The remaining variable u 2 will enter, (along with u 2 , u 2 ) explicitly in the system and is 
therefore advisable (if not mandatory) to be used as the time parameter, i.e. to effect 
the change of time coordinate 

t -> u 2 {t) = s, ui(t) -> Ui(t(s)), u 4 {t) -> u 4 (t(s)), ue(t) -> u e (t(s)). (3.27) 

This choice of time will of course be valid only if u 2 is not a constant. We are thus 
led to consider two cases according to the constancy or non-constancy of this dependent 
variable. 

Until now, we haven't commented upon the range of values that the parameter h 
can attain. As it is well known, for the value h = we come across the Class A model, 
which admits a Lagrangian description, whereas for h 7^ we have the Class B model 
which lacks such a Lagrangian description. So we are forced to examine two further 
possibilities, as to whether h is equal to, or different from, zero. 

3.1.1 Case I: h — and u 2 (t) = k 2 

In the parametrization (13.25}) the determinant of y a p, is 

^W4 e3 " 1+4U6 (l-4fc 2 2 ) 
so we must have — | < k 2 < | . The third linear constraint reads 

8 k 2 Ua 

E 3 = 0^ 2 = u 4 = k 4 or k 2 = (3.28) 
— 1 — I — 4 k 2 

The case u 4 = k 4 leads, through equation E u = to k 2 = . Thus, the only possibility 
is k 2 = 0. Substituting this value into the quadratic constraint equation E we obtain 

-- (3m 1 2 + 8m 1 m 6 + 4m 6 2 ) = (3.29) 
which has the following two solutions 

u\ = k\ — 2uq (3.30a) 
2 

ui = ki uq (3.30b) 

3 



s 



For the first of (I3.30p all the spatial EFE's are equivalent to the equation 

2 uq NN + 2 N 2 (uq — ii 6 ) = (3.31) 
from which we have for the lapse function 

N 2 = ke -2u a{i 2 (33 2 ) 

Choosing a time parametrization uq = — | ln(^), and using the automorphism ma- 
trix (13.71) with entries Si = \ (In 2 — k\) , s 2 = S3 = S4 = we arrive at the line element 

ds 2 = -dr 2 + r 2 dx 2 + dy 2 + dz 2 (3.33) 

which describes a flat space admitting a manifest VIIq symmetry [13J. To the best of our 
knowledge, it is the first time that this line element emerges in the course of investigation 
of the solution space to this Bianchi Type. 

For the second of (I3.30p all the spatial EFE's are equivalent to the equation 

2ii 6 NN -2N 2 (u 2 6 + u 6 ) =0 (3.34) 

which gives the lapse function 

N 2 = ke 2u6 iil (3.35) 

Choosing a time parametrization uq = | In (4^), redefining the constant k\ = | In |^ 
and using the automorphism matrix (13.71) with entries s x = | In s 2 = s 3 = s 4 = we 
arrive at the line element 

ds 2 = -t dr 2 + - dx 2 + t 2 dy 2 + t 2 d z 2 (3.36) 

T 

This line element was first derived by Ellis [14j.[15] and admits, besides the three 
killing fields ( 13.31) (with h = 0), a fourth symmetry generator 

£4 = d x (3.37) 

along with a homothetic vector field 

r] = 2rd T + 4:xd x + ydy + zd z (3.38) 

There is thus a G4 symmetry group acting (of course, multiply transitively) on each V3 
of this metric. However, it is interesting to note that we have not imposed the extra 
symmetry from the beginning, but rather it emerged as a result of the investigation 
process. 
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3.1.2 Case II: h = and u 2 (t) = t 

With this choice of time gauge the third linear constraint reads 



E* = 



it* 



U4 = =r- W4 = /C4 



4t 2 - 1 

With this information at hand, the quadratic constraint _E yields the lapse function 



(3.39) 



U\ 



N 



64 1 2 



((4 1 2 - 1) (ux + 2 rifl) (3 ux + 2 ri 6 ) + 16 1 {u x + ri 6 ) + 4) 



(3.40) 



We now turn to the spatial equations of motion and substitute the above lapse. The 
simplest is E 33 = and the coefficient of ii x in this equation is proportional to the 
quantity 

((4t 2 - l)u 6 + 2t) Mi + (4t 2 - I) ill + 4tw 6 + 1 

which can be safely regarded different from zero, since by setting this quantity equal 
to zero and solving for iii we end up with zero lapse (with the help of the rest of the 
equations of motion). We can thus solve -E33 = for u\ and substitute into En = 0. In 
this transformed equation En = 0, the coefficient of uq is proportional to 

(-2 1 + sin(2 fc 4 )) (tti + ttfi) - 1 

a quantity which is different from zero, since it's nihilism leads again to zero lapse. From 
the transformed En = we have the expression for ii 6 , so we finally arrive a the following 
polynomial system of first order in ui, m 6 



ui = (ui\ Ax \u 6 ), ii 6 = (ux\ A 2 \u 6 ) 



(3.41) 



where we have introduced the notation (ttj 
matrices Ai, A 2 given by 



luiufuf) and \ui) = with the 4x4 



/ 



4t 2 -l 
4t 2 +l 
t(4 t 2 -l) 



3(l-4t 2 ) 
4t 



16 1 
4t 2 -l 

4 

2(-4f 2 +l) 
t 





4 0\ 
^1 








0/ 



.4, 



-28 t 2 +l 
4t 2 +l t(4i 2 -l) 
16 f _i 9 

-4t 2 +l ±Z 
o -12f 2 +3 



jt 2 +2 



V 



If 




f 






-4f 2 +l \ 

t ' 





y 



(3.42) 



Due to the form of An A 2 (their components are rational functions of the time t), system 
(I3.4ip can be partially integrated with the help of the following Lie-Backlund transfor- 
mation 



us(t) 



-16r 2 (t) + 1 



-Atr(t) 



4(4t 2 - l)r(t) 

-16r 2 (t) + 16tr(t) + 3 
8(-4t 2 + l)r(t) 
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(3.43) 



+ 2tf(t) 



resulting in the single, second order ODE for the function r(i) 

1 

r 



„ (12t 2 + l)r + t 
•- r 2 + —, — — — '—. — r 



At 2 + l)tr 

At this stage, in order to solve (13.441) we apply the contact transformation: 

CMO , i i ( w (0-i)(^(0 + i) 3 



r(t) = - — , t = - 
4 MO 2 



which reduces it to 



w"(0 



1 + 10(0' r(t) s™'(0M0 
«/(0 2 <M0 i MO 2 - 1 



(3.44) 



(3.45) 



(3.46) 



MO e 2 f 

which is nothing else but the third Painleve transcendent w := Pm(a, /3, j, 5) with 
entries (a, (3, 7, 5) = (—1/2, 1/2,0,0). For completeness we give the general form of the 
equation that the third Painleve transcendent satisfies: 

*® = lM--t + -( + ' ),al€ > + «i5 (3AT> 

Using the final equation (13.461) . the contact transformation (13.451) and the Lie-Backhand 
transformation (13.431) we find that the functions Uq and u[ are given by 



MO 



In 



MO 



'MO + 1) 2 
1 



MO 



MO = t^% + \ w (0 + A „. , 

4w{£,) z 4 4w(£) 44 2 



and the lapse function has the form 

N 2 = 

The scale-factor matrix 7 Q/ g is thus 

(MO + i) 
(MO - 1) 



lap 









MO 










V 






iee 



3 "1 



(3.48a) 
(3.48b) 

(3.49) 



MO 



MO 



(MO - 1) 
(MO + 1) 












(3.50) 



which can be brought to diagonal form with the aid of the automorphism matrix ( 13.71) 
with entries S\ = —\ In 2, s 2 = —^5 



la(3 



MO 






MO 







i 

w(0 










\ 



(3.51) 
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Gathering all the pieces we arrive at the final form of the line element 

1 



ds' 



16 4 4 



/ill < 



w sin 2 x + cos 2 x) d y A 



+ 



sin(2 x) (w — 1) 



dydz+\j — (w cos 2 x + sin 2 x) dz 2 ^ (3.52) 



which represents the general solution of Bianchi Type VIIq non-flat Vacuum Cosmology, 
since it contains the expected number of three essential constants (two implicit in the 
third Painleve transcendent plus the overall k). The above line element was first given 
by Lorenz-Petzold [16], but it was not then pointed out that it represented the general 
solution. 



Particular Solutions 

In order for the contact transformation (13.451) to be well defined, it is obvious that 
the function must not be constant. However, remarkably enough, the resulting line 
element (13.521) does not inherit this restriction. Thus, if there is some constant solution 
to equation ( 13. 46ft . it could produce a particular solution through ( 13.521) . By inspection 
it is obvious that (13.461) admits the solutions w(£) = ±1, so we could use them to obtain 
two particular solutions. 

► Subcase w(£) = 1 

With this value of u>(£) ( 13.521) indicates that £ > 0, so plugging this value to ( 13.481) 
we have 

ui(0 = *i-i]ii£ we(0 = -±fci + ±ln(16 £ 3 ) (3.53) 

which, after using the usual simplifications brought by the automorphism matrix (13.71) 
and redefining the variable £ to £ = r 4 , results in 

ds 2 = -rdr 2 + -dx 2 + r 2 dy 2 + r 2 dz 2 (3.54) 

T 

which is the line element ( 13. 36ft . 

► Subcase w(£) = —1 

Now from (13.521) we must have £ < and from (13.48bl) we obtain 

ui(0 = *4-£-|H£| (3.55) 

while from (I3.48al) u 6 remains undefined. The line element (13.521) with the help of the 
automorphism matrix (13. 7p and the definition £ = — r 4 becomes 

ds 2 = k 2 (e r4 rdr 2 H e 1 " 4 dx 2 + r 2 cos(2 x) dy 2 — 

V 4r 

2 r 2 sin(2 x) dydz — t 2 cos(2 x) dz 2 ) (3.56) 
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which, though physically acceptable, corresponds to Bianchi Type VIIq symmetry on 
T 3 , and was first given by Barnes [TT] . 



3.1.3 Case III: /i / and u 2 (t) = k 2 

In this case the determinant of the scale factor matrix 7 a/ 3 is 

det( 7af3 ) = ~ e 3ui+4 ^ (1-4 k\) (3.57) 

so we must have — \ < k 2 < \ in order for it to be positively defined. 
The third linear constraint £"3 = reads 

1 4 A; 2 

4fc 2_ I {8k 2 2 u 4 + 2h(l-4k 2 2 )u 6 ) =0^u 6 = k 6 + h ^ k 2_ X) u ± ( 3 - 58 ) 

and the quadratic constraint E a = gives for the lapse function N 2 

N " = 16 h* (4A; 2 -l)(3/i 2 (4k 2 -l)-4k 2 ) ( 3 h " (4 k2 ~ 1)2 ^ + 

32 h k\ (4 k\ - 1) u x ii A + 16 k\ (h 2 (4 k\ - 1) + 4 A; 2 ) ?i 2 ) (3.59) 



Now we are ready to attack the spatial equations of motion after substituting in them 
the above lapse. -E33 = is again the simplest one. In this equation, the coefficient of 
U4 is proportional to 

k\ mi (h (4 k\ - 1) u x + (4 k\ + h 2 (4 k\ - 1)) u 4 ) (3.60) 

so in order to solve _E 33 = for ii 4 we must ensure that the above quantity is different 
from zero. Setting this quantity equal to zero we get 

Ak 2 + h 2 (4fc 2 -l) 
ux = M4fei -i) ^ ( 3 ' 61a ) 

Ui = fcj (3.61b) 

A; 2 = (3.61c) 

► The solution ( 13. 61 aft leads to inconsistency. 

► The solution (13.61bjl forces equation E 33 = to give either k 2 = ±-^—^=== which 

leads to zero lapse or ko = ± , which makes the determinant of 7„« negative so 

P 2 V / ^1 
is unacceptable. 

► The solution (13.6 left satisfies all the spatial equations and leads, after the usual 
simplifications achieved by the automorphism matrix ( 13.71) and the choice of time gauge 
U\ = 2 ln(2 hr), to the line element 

d s 2 = _ d T 2 + h 2 T 2 dx 2 + e 2h , T 2 dy 2 + e 2h , ^ rf ^ ( g ^ 
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which describes a flat space admitting a manifest VI Ih symmetry, a line-element first 
presented by Doroshkevich et al [IB] and reproduced by Siklos [T§] . 

Having ensured that the term (13.601) is not equal to zero we can solve E 33 = for ii 4 
and substitute into the other equations of motion. From En = we have 

2{3h 2 (4 k\ - 1) + Ak 2 ) 



(sin(2 u 4 ) -2k 2 )(u l + 3^31) 



WW^T) "'' =0 (3 - 63) 



which leads to the following possibilities 



U4 = - arcsin(2 k 2 ) (3.64a) 



- = 3MVl) "« <3 " 4b) 
2(4/L 2 (/i 2 -l)-/i 2 ) , , . 

ui = fci+ l) AU2 — - -u 4 (3.64c) 



2(3/i 2 (4/c| - 1) +Akf 

3/i (4/c| — 1) 
2(4 fc 2 (fr 2 - 1) - fr 2 ) 
M^f-l) 

► The solution (I3.64al) leads to zero lapse. 

► The solution ( 13.64bl) leads to inconsistency. 

► The solution (I3.64cl) satisfies all the spatial equations and leads, after the usual 
simplifications with the automorphism matrix (13.71) and the choice of time gauge Ui = r, 
to the line element 

"1 2 /i^ (-^^ 1) 

ds 2 = -e T (-dr 2 + dx 2 ) -2\e 2h{T+x) sin2(r + x) dydz+ 

e 2h{T+x) ^1 + A cos2(r + x)^ dy 2 + e 2h{T+x) (l - A cos 2(r + x)J dz 2 (3.65) 

which was presented in [18] and [T9]. This line element admits, besides the three killing 
fields (13.31) . three more, namely 

-A 2 + fe 2 (A 2 -l) , _ s -A 2 + fe 2 (A 2 -l) , _ x 

£ 4 = e ma 2 -d T >d T -e ma^-d (XT) d x (3.66a) 

-A 2 + fe 2 (A 2 -l) f , -A 2 + ft 2 (A 2 -l) , , 

£ 5 = ye ma 2 -d I s r ^ r - ye ma 2 -d ^ r ^ 

^2 ^2 ^"j 

_l_ e ma 2 -i) (x+r) ( Cl C os 2 (r + x) + c 2 sin 2 (r + x) + c 3 ) d y 

^2 ^2 \\ 

+ e ma 2 -d (x+r) ( C2 C os 2 (r + x) - ci sin 2 (r + x)) <9 Z (3.66b) 



-A 2 + ft 2 (A 2 -l) , _ \ -A 2 + ft 2 (A 2 -l) , _ \ 

ze ma 2 -d [x T >d T -ze ma 2 -d (x T *d x 

+ e ma 2 -d (x+r) ( C2 cos 2 (r + z) - ci sin 2 (r + x)) 

-A 2 -fc 2 (A 2 -l) , v 

_ e ma 2 -d 1 + j (ci cos 2 (r + x) + c 2 sin 2 (r + x) - c 3 ) <9 2 (3.66c) 



14 



where the constants (ci, c 2 , c 3 ) are given by 

Xh(X 2 + h 2 (X 2 -l)) 
4(A 4 + /i 4 (A 2 -l) 2 + 2/i 2 (A 2 -l)(3A 2 -2)) 1 ' ; 

Afr 2 (A 2 -l) 

2 2 (A 4 + /i 4 (A 2 - l) 2 + 2h 2 (X 2 - l)(3A 2 -2)) l ' J 

C3 = 4 (A 2 + Aa 2 - 1)) (3 - 6?C) 

Again it is worth mentioning that this Gq symmetry was not imposed from the begging 
but emerged during the seeking of the solution space. The non-vanishing commutators 
are 

[&,6] = 6i = & 

Ks, &] = 2{^ + h)U fc] = 2 (SL + £) £ 5 + £ 6 

C2 C2 / 

[6,ee] = -6 + 2(- + J)e6 (3.68) 

with (01,02) given by (13.671) . Finally the line element (13.65!) admits a homothetic vector 
field 

71 = -A 2 +W-l) ^ - -A 2 + W- 1) O'+V 9 *** 9 * ( 3 - 69 ) 
3.1.4 Case IV: h ^ and u 2 (t) = t 

In this case the determinant of the sale factor matrix becomes 

det[ la p] = - A e^ +4u «(l-At 2 ) 

so we must demand that |i| < - in order for to be positive definite. 
The third linear constraint E3 = can be used to define the function uq 

2t(2tu 4 -h) , [2t(2tu 4 -h) , . „ . 

The quadratic constraint £7 = defines the lapse function iV 2 

^ = 16tf(4 < »-l)(3W(4t»-l)-4 < ') ( 3 " 2 (4 ' 2 " 1)2 * ! + 

32 /i t 2 (4 1 2 - 1) m 4 + 16 £ 2 (/i 2 (4 1 2 - 1) + 4 1 2 ) m 2 - 4 /i 2 ) (3.71) 
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Substituting the above values of the lapse N 2 and the function w 6 in equation E 33 = 
we find the coefficient of U\ to be proportional to 

At 2 {At 2 + h 2 (At 2 - 1)) ill + Aht 2 (At 2 - I) U!U 4 - h 2 

a quantity that can be safely regarded different from zero, since it's nihilism leads either 
to zero lapse or to inconsistency. Thus we can solve £33 = for iii and substitute it to 
En = 0. In order to solve this equation for it 4 we must be assured that it's coefficient 
does not vanish. Setting this coefficient equal to zero we arrive at the following equation 



Ui 



hcos2u 4 + 2t (At 2 (h 2 + 1) - h 2 )u 4 - 2ht 



2t(sm2u 4 -2t) 



h(At 2 -V 



which is unacceptable because it leads to inconsistency. After solving equation En 
for M4 we finally arrive to the following polynomial system of first order in iii, it^ 



Hi = (ui \ Bi \u&), U4 = (ui \ B 2 1^4) 



(3.72) 



where we have used again the notation (iii 
matrices Bi, B2 given by 



luiufuf) and {iii) = (ui\ with the 4x4 



Bo 



/ 4(h 2 -4(-l+/i 2 )t 2 ) 
: 



Bi =f 







(l-4t 2 ) 4 

-4t(l-8f 2 +6fc 2 (-1+4 t 2 )) 32 t 2 (-fe 2 +4(-l+fe 2 )f 2 ) 
(-1+4 1 2 ) 3 h(-l+4 t 2 f 

32 t 3 



-3fe 2 +12(-l+fc 2 )t 2 

(l-4t 2 ) 2 
3t 



-1+4 1 2 



2h 



4t 2 +h 2 (l-4t 2 ) 



3fe(l-4t 2 ) 2 
-8t 2 +h 2 (-2+8 1 2 ) 





h{i-it 2 Y 




-6h 2 (l-4f 2 ) 2 +4f 2 (-3+8 1 2 ) 
-it 3 +h 2 t(-l+it 2 ) 
16 t 2 (-1+4 1 2 ) 
-h 2 +A(-l+h 2 )t 2 
3t(l-4f 2 ) 2 



-h 2 +i(-l+h 2 )t 2 





16f 2 (-16t 4 +fc 4 (l-4t 2 ) 2 ) 
h? (l-4i 2 ) 4 
-At (-fe 2 +4(-l+fe 2 )f 2 ) (4 t 2 +3h 2 (-1+4 1 2 )) 
h 2 (-1+4 t 2 ) 3 









o\ 




0/ 



;(4 t 4 +h 2 1 2 (-1+4 1 2 )) 

-4ht 2 +h 3 (-1+4 1 2 ) 

32 t 3 (-1+4 1 2 ) 
~4ht 2 +h 3 (-1+4 1 2 ) 






-4t(-3fe 2 +4t 2 +12fe 2 t 2 ) \ 
h 2 ' 








where 



./ 



3h 2 (At 2 



(At 2 - l) 2 



At 2 (At 2 - I) (3 h 2 (At 2 -I) -At 2 
, 9 = : 



A(h 2 -l)t 2 - h 2 



(3.73) 



Due to the form of Bi, B 2 (their components are rational functions of the time t), system 
(I3.72p can be partially integrated with the help of the following Lie-Backlund transfor- 
mation 



Ui 



u.i 



h 2 (3h 2 (At 2 - 1) +At 2 ) tanr(t) 4 (h 2 - 1) t 2 - h 2 . 



2hty/l-lt 2 (3h 2 (At 2 - 1)-At 2 ) 2h^l-At 2 



3/iVl -4t 2 tanr(t) 



At(3h 2 (At 2 



4t2) + i v / T^r(t) 



r(t) (3.74a) 



(3.74b) 
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yielding the single second order ODE for the function r(t) 

^ Vl -At 2 + tanr^ r 2 + 

3 (4 t 2 (3 h 2 - 1) - h 2 ) h (3 h 2 (4 1 2 - 1) - 8 t 2 ) 

tanr f + 



t (4 (3 h 2 - 1) t 2 - 3 /i 2 ) t (4 (3 h 2 - 1) t 2 - 3 /i 2 ) VI -4t 2 

9/i 4 (sin2r + /i v / l-4t 2 ) 2 /i (3 /i 2 (l - At 2 ) + 8t 2 ) 2 

sec r 



2t 2 (At 2 + 3 /i 2 (1 - 4t 2 )) 2 2t 2 (3/i 2 (1 - 4t 2 ) + 4t 2 ) 2 ^(1 -4t 2 ) 3 

(3.75) 

This equation contains all the information concerning the unknown part of the solution 
space of the Type Vllh vacuum Cosmology (h ^ 0). Unfortunately, it does not posses 
any Lie-point symmetries that can be used to reduce its order and ultimately solve it. 
However, its form can be substantially simplified through the use of new dependent and 

independent variable (p, w(p)) according to r(s) = iarcsin— ^==, s = yj 12 3 fc 2 

thereby obtaining the equation 

h(l-il 2 ) -2 h 2 (l-il 2 ) 2 
U = ± — , => U = - — — ; : . H3.76) 

^(<oh 2 p + A-Qh 2 ) (p 2 -u 2 -l) (6h 2 p + A-6h 2 ) (p 2 -u 2 -lf 1 

This equation is a special case of the general equation 

u 2 - ^-^ 2 (3 77) 

{K + \p){p 2 -U 2 -l) 1 } 

4 

with the values k = — 6 + — , A = 6. The general solution of (13.771) was first given in 

h l 

and can be obtained as follows: First we apply the contact transformation: 



( \ 8 . 4 ( 2 ^~ 1 ) >(t\ k A 

Hp) = -jy{0 + ^ — y{0 p = ~\ + \ y w 

w(p) = 2e-i hp) 



(3.78) 



2y"(0 

which reduces it to 

e (£ - l) 2 2/" 2 = - V (£ y' - y) 2 + 4 y' 2 (£ y' - y) - ~ y' 2 + y' (3.79) 

This equation is a special form of the equation SD-Ia, appearing in [2T] . where a 
classification of second order second degree ordinary differential equations was performed. 
The general solution of (13.791) is obtained with the help of the sixth Painleve transcendent 
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w := Pvi(ct, P, 7, 5) and reads: 

e 2 (C- I) 2 



4w(w-l)(w-^) 



w — 



w (w — 1) 



+ -(l±\/2^) 2 (l-2w)- 



7 
4 



1 - 



2(£-l) 



4 

1 5 



2£ 



2£(™-l) 



where the sixth Painleve transcendent w := Pvi(ct, /3, 7, 5) is defined by the ODE: 



(3.80) 



w 



1 



2 \w — 1 u> If — £ 
w (w — 1) (w — £) 



./2 



1 1 
- + 



e (e-i) 5 



a + + 7 



+ 5 



if!- 



(3.81) 



The values of the parameters (a, P, 7, 5) of the Painleve transcendent, can be obtained 
from the solution of the following system: 

a-/3 + 7-5±v / 2^+l = ~ (3.82a) 



09 + 7) (a + tfiv^) = (3.82b) 



k 2 - A : 



(- - V) ( n - r) ± v 2d : I ) + 1 (a - - 7 + 5 ± v 7 ^) 2 ]G 



- (7 - /?) (a + 5 ± v 7 ^) 2 + 1 (/3 + 7) 2 (a - S ± V2^ + l) 







(3.82c) 
(3.82d) 



Plugging in (13.821) the values of k = — 6 + — , A = 6 for Type Vllh, we have 

hr 

twenty-four solutions (counting multiplicities) of this system. In order for the parameters 
(a, P, 7, 5) to be real numbers we end up only with four possibilities 



(a,P,~f,6) 

(a,P,-f,5) 

and 

(a,/3,7,<y) 

(a,/3,7,5) 



Ah 2 - 1 
2/i 2 

4/i 2 - 1 
2/i 2 



1 1 



1 l-2h 2 



+ \ 3 



/i 2 '2/i 2 ' 2/i 2 ' 2/i 2 
1 1 - 2 /i 2 



1 1 



/i 2 '2^ 2 ' 2/i 2 ' 2/i 2 



v/3 



(3.83a) 



1 2-3h 2 Vl-3/i 2 3/i 2 -2 VI - 3 /i 2 1 
2' 2/i 2 + ^ ' 2h 2 + *2 '2 



1 2-3/i 2 _ VI -3/i 2 3/i 2 -2 Vl-3/i 2 1\ 1 



(3.83b) 
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For the values h = ±— = the above relations coincide and as we will show these values 

of h give rise to a particular solution. 

Gathering all the pieces the final form of the general line element describing Bianchi 
Type Vllh vacuum Cosmology is 



2 ,os(2» j (u)\/e(e-i)(y / (0-^)^ 1 ^ 2 



where 



+VW^) (Vm+sm(2u,m - ^\ (<r 2 ) 2 + e-(« (<r 3 ) 2 ) (3.84) 



, m (-1 + /1 2 ) (-l + 2Q + 2fe 2 y(0 



M4(0 = 4 M -i + 0e(-i + ^^)) (3 - 85b) 

and ?/(£) is given by (13.801) . Again, this line element contains three essential constants, 
thus representing the general solution to the EFE's for the Class B Vllh case. 

Particular Solutions 

Even though the line element (13.841) represents the general solution of Bianchi Type 
Vllh vacuum Cosmology, it does not come into a manageable form due to the appearance 
of the sixth Painleve transcendent. To partially remedy this inconvenience, we give, in 
the following, some closed form line-elements arising from particular solutions to (13.801) 
and fl3~%Tj) . 

► Subcase y(£) = c and \h\ < — -= 

One way to obtain a particular solution from the above line element (13.841) is to 
follow the reasoning of Case II, i.e to observe that, although the form of the contact 
transformation (13.781) implies that the function y(£) cannot be constant, the line element 
(13.841) is free of this restriction; the difficulty with the negative argument in the square 
root is circumvented by using the hyperbolic sine/cosine (see (13.861) below). We can thus 
check if the assumption ?/(£) = c leads to a particular solution. Skipping the calculational 
details, we find that for y{£) = c = ^ 1 2 ~ h 3 2 h ' 2 all the Einstein's field equations are satisfied 
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and we end up with the new line element 

ds 2 = k 2 sin(4/ir)(/(r)(dr) 2 + sin(/ i In /(r)) (cr 1 ) 2 - sin(/i In /(r)) (cr 2 ) 2 

+2 cos(/i In /(r)) cr 1 cr 2 + /(r) (cr 3 ) 2 ) (3.86) 



/(r) = sin~^(4/ir) tan _y ^5^(2/ir), < -7= 

v3 

which even though is physically acceptable it corresponds to Bianchi Type Vllh sym- 
metry on T 3 . Since the above line element admits only the three killing fields (13.31) and 
no homothetic vector field we can conclude that the constant k is essential. ^ 
An interesting property of the line element ( 13.861) is that, for the value h 2 = -, i.e. 
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as = k (esc — (ar) — sin — p sin (v 3 In sin — p J (cr ) 
^ v3 v3 v3 

4 T / j 4 T \ , o \ o 

+ sin —= sin (V3 msin —7= J (cr ) 
V 3 v 3 

+2 sin ^ cos (V3 In sin -^L) cr 1 cr 2 + esc 2 -^(cr 3 ) 2 ) (3.87) 
v3 v3 v3 ' 

admits a fourth killing field, namely 

-2x 4 T -2 g 4 T 

77 = e vs s in — 9 T -2e^3 cos — p d x (3.88) 
v 3 v 3 

The geometry (13.871) was first given by Petrov [22] and it is the only vacuum solution 
admitting a simply transitive G4 as its maximal group of motions. This group of motions 
has two subgroups G3 of Bianchi Types I and VII h 2 = i acting in time-like hyper-surfaces. 

► Elementary solution of Painleve transcendent 

As it is well known, although for generic values of the parameters (a, (3, 7, 5) the 
Painleve functions are transcendental, there exist a lot of elementary solutions for spe- 
cial values of these parameters [25] . [24] . In the case at hand the following Lemma is 
applicable 

Lemma The function w satisfying u>(£) 2 — 2£u>(£) + £ = is a solution of (13.811) 
when the parameters (a, /3, 7, 5) obey the relations a + 8 — ^, (3 — —7. 

Proof Direct computation. □ 



2_ 

from the first of (13.831) we have (a, /3, 7, 5) = (~, y , — y , |). Choosing now the parametriza- 



Using (I3.83p . the conditions of the above Lemma are fulfilled for h = i~7fj- Then 



1 11 _n 3 

>8' 8 ' 8 ' 8 

tion 



w(0 = ^~ 2 e 4hT , £ = cosh 2 (2 /it) (3.89) 
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we can compute y(r) from (I3.80P (with the minus sign) and m^t), -u 4 (t) from (I3.85p . 
thereby arriving at the following line element 

ds 2 = k 2 ^-eir sinh-i(4/ir)(dr) 2 + e- 2/lT sinh^(4/ir) (e 4/lT + sin(4r)) (cr 1 ) 2 

+2 e- 2hT cos(4 r) cr 1 <x 2 + e~ 2/lT sinh^ (4 h r) (e 4?iT - sin(4 r)) (<r 2 ) 2 

sinh-i (4 h r) (cr 3 ) 2 ) (3.90) 

This geometry was first given by Lukash [23] and, like (13.861) . admits only the three 
killing fields (13.31) and no homothetic vector field. Therefore, the constant k is essential. 



4 Discussion 

In [7], [20] a systematic approach for investigating the solution space of Bianchi Type 
Cosmologies was developed by the use of automorphisms and the theory of symmetries 
of ordinary, coupled differential equations. The result was the comprehensive recovery 
of all known closed form Type III solutions, as well as the presentation of the general 
solution in terms of the sixth Painleve transcendent. In the present work we have applied 
the method to the case of Bianchi Type Vllh family of vacuum geometries. Again, the 
general solution is implicitly given in terms of the third (I3.52j) Painleve transcendent 
or the sixth Painleve transcendent (13.841) for the Class A (h = 0)) and the Class B 
(h ^ 0) case respectively. Through the investigation of either Particular or Elementary 
solutions of the Painleve transcendents we are able to concisely recover, in a systematic 
fashion, all six known solutions (l3^ . f[33oT) . (l3^ . f[3^ . fl3^ All these metrics 

have originally been obtained in a time scale of 20 years or so, by prior assumption of 
symmetry and/or other physical requirements; e.g. Petrov' s solution [22] was derived 
with the use of automorphisms seeking G4 homogeneous metrics while Lukash' s solution 
[23] was derived based on a physical interpretation of Type Vllh cosmological models, in 
terms of circularly polarized gravitational waves of arbitrary wavelength in a space having 
constant negative curvature. Their reacquisition single-handed, proves, we believe, the 
value of our method. A very important result is, of course, the discovery of the new 

family of solutions (I3.86P for the range of the group parameter h 2 < -. Besides of the 

3 

obvious value of a new family of solutions to the EFE's it also points to the unexpected 
existence of a sector with particular behavior for this Bianchi Type. It is known that Type 

Vlh model has an exceptional sector corresponding to the value h 2 = - but, for Type 

9 

Vllh such a behavior is first observed. The fact may be taken as a further strengthening 
evidence of the widespread belief that the two Types are very much similar. We hope 
that the application of the method to Type Vlh will bear analogous fruits. As for Types 
VIII, IX, the recent discovery that some particular configurations are described by the 
third Painleve transcendent [26] strengthens our belief that their solution space will also 
be attained by our method. We plan to return to these issues in the immediate future. 
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Finally, we deem it useful to end this discussion by briefly describing the investigated 
solution space through the following tables: 
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Bianchi Type VIIp metrics 



Line Element Isometry Type Comments 

ds 2 = -(dt) 2 + ((T 1 ) 2 + (a 2 ) 2 + r 2 (a 3 ) 2 G 10 onV 4 Flat Space 



G 4 on V 3 , t > 



(is 2 = -r(dr) 2 + r 2 ( C r 1 ) 2 + r 2 ( - 2 ) 2 + -( C r 3 ) 2 G 4 on T 3 , r < LRS 

r 

ds 2 = K 2 ( e r4 r (dr) 2 - r 2 (cr 1 ) 2 + r 2 (<t 2 ) 2 + — (<r 3 ) 2 ) G 3 on T 3 , r > Non-homothetic 



MO 



16f 



rf s 2 = K 2 _<L^ ( d ^2 + (o- 1 ) 2 + i / ( CT 2 ) 2 + e " l(f) (<x 3 ) 2 G 3 on V 3 General Solution 



™(0 



where the 1-forms cr a are given by 



er 1 = sinxdy + cosxdz, cr 2 = cos xdy — sin xdz, er 3 = -dx (4-1) 



and is defined by equation (I3.48bl) 

4w(0 2 4 4w(0 4£ 2 

with standing for the third Painleve transcendent w := Pm(-i ~, 0, 0), defined by fl3~4T3j) 



Bianchi Type VHh metrics 



Line Element 


Isometry Type 


Comments 


ds 2 = _( dr )2 + h 2 T 2 (<T 1)2 + h 2 T 2 ^2)2 + 4 ^2 r 2 ^3)2 


G\o on V4 


Flat Space 


ds 2 = \ exp (- 2A2 + ( y 2 g 2_1) r) ( dr 2 + 4(<x 3 ) 2 ) + |e 2 ^ (1 + A sin2r) (<x 2 ) 2 
-ie 2ftT (-1 + A sin2r) (cr 1 ) 2 + e 2ftT A cos 2r (cr 1 ) 2 (cr 2 ) 2 


G 6 on V A 


Homothetic 


rfs 2 = k 2 sin(4 hr)[f{r) (dr) 2 + sm(h ln/(r)) (cr 1 ) 2 - sin(/i ln/(r)) (cr 2 ) 2 
+2 cos(/i In /(r)) cr 1 cr 2 + /( r ) (cr 3 ) 2 ) 


G 3 on T 3 


Non-Homothetic 


rfs 2 = /t 2 ( esc 2 ^§(dr) 2 - sin ^| sin (y/3 In sin ^) (cr 1 ) 2 
+ sin % sin (V3 In sin %\ (cr 2 ) 2 
+2 sin ^ cos (v/3 In sin cr 1 cr 2 + esc 2 ^(o" 3 ) 2 ) 


G 4 on T 3 


Non-Homothetic 
Maximal G4 

h 2 = \ 


cis 2 = k 2 ( - e^r sinh"l(4/ir) (dr) 2 + e~ 2hr sinh^(4/ir) (e 4fer + sin(4r)) (cr 1 ) 2 
+2e~ 2ftr cos(4r) cr 1 cr 2 + e' 2hr sinh5(4/ir) (e 4hT - sin(4 r)) (cr 2 ) 2 
+e^r sinh"i(4/ir) (cr 3 ) 2 ) 


G 3 on V3 


Non-Homothetic 
h 2 - ± 


ds 2 = k 2 ^ - -^0^ (d£) 2 + y/Z (£ - 1) (vMO - sin(2 ^(0 " 

+2 cos(2 M4 (0) (£ - 1) (y'(0 - ct1 ct2 

+ (e - 1) (vV(o + «m(2 ^(0) vV(o - £) (^ 2 ) 2 + e " l(f) ( ct3 ) 2 ) 


G 3 on V3 


General Solution 



where the 1-forms <r a are given by 
^.l _ e hx ( s { nx( iy _)_ cosx dz) , er 2 = e hx (cos a; dy — sin a; dz) 
the function /(r) stands for 

1 Vl-3 

/(t) = sin~^ (4 /i r) tan ^~ (2 h r) 
and the functions Wi(£)> M 4(0 are defined by (13.851) 

(-1 + /1 2 ) (-l + 20 + 2/i 2 i/(0 



<(£) 

with y(0 defined by (|£B0|1 . 



l-2£ + 2/i 2 y(0 



4/i(-l + £H(-l + />V(0) 



25 



References 

[1] O.Heckman and E. Schiicking, Relativistic Cosmology in Gravitation (an introduc- 
tion to current research) edited by L. Witten, Wiley (1962) 

[2] A. Harvey, Jour. Math. Phys. 20, 251 (1979) 

[3] T. Christodoulakis, G. O. Papadopoulos and A. Dimakis, J. Phys. A 36, 427 (2003) 

[4] R.T. Jantzen Comm. Math. Phys. 64, 211 (1979); JMP 23, 1137 (1982); C. Uggla, 
R.T. Jantzen and Rosquist, Phys. Rev. D 51, 5522 (1995) 

[5] J. Samuel and A. Ashtekar, Class. Quan. Grav. 8, 2191 (1991) 

[6] T. Christodoulakis, G. Kofinas, E. Korfiatis, G.O. Papadopoulos and A. Paschos, 
JMP 42, 3580 (2001) 

[7] T. Christodoulakis and Petros A. Terzis, J. Math. Phys. 47, 102502 (2006) 

[8] G.F.R. Ellis and M.A.H. MacCallum, Commun. Math. Phys. 12, 108 (1969) 

[9] "Exact Solutions of Einstein's Field Equations" (Second Edition), H. Stephani, D. 
Kramer, M. MacCallum, C. Hoenselaers and E. Hertl, Cambriodge Monographs on 
Mathematical Physics, CUP, Cambridge (2003) 

[10] T. Christodoulakis, E. Korfiatis and G.O. Papadopoulos, CMP 226, 377 (2002) 

[11] "Differential equations: Their Solutions using Symmetries", H. Stephani, Edited by 
M.A.H. MacCallum, Cambridge University Press, Cambridge (1989) 

[12] See e.g. "Applications of Lie Groups to Differential Equations", Peter J. Olver, 
Graduate Texts in Mathematics 107, (2000) 

[13] A.H. Taub, Annals of Mathematics 53, 472 (1951) 

[14] G.F.R. Ellis, JMP 8, 1171 (1967) 

[15] J.M. Stewart and G.F.R. Ellis, JMP 9, 1072 (1968) 

[16] Lorenz-Petzold, Acta Phys. Polon. B 15, 117 (1984) 

[17] A. Barnes, J. Phys. A 11, 1303 (1978) 

[18] Doroshkevich A.G., Lukash V.N. and Novikov I.D. Sov. Phys. JETP 37, 739 (1973) 

[19] S.T.C. Siklos, J. Phys. A: Math. Gen. 14, 395-409 (1981) 

[20] T. Christodoulakis and Petros A. Terzis, CQG 24, 875 (2007) 

[21] Christopher M. Cosgrove and George Scoufis, Stud. Appl. Math. 88:25-87 (1993) 

26 



Petrov A.Z., Gravitational field geometry as the geometry of automorphisms in 
Recenet developments in general relativity, Pergamon Press-PWN Warsaw, Oxford 
(1962) 

"The Analytic Solutions of the Painleve Equations", Gromak V. A. and Lukashevich 
N. A, Minsk: Universitetskoye Publishers (1990) (in Russian) 

[24] http:/ /www.newton.cab.ac.uk/webseminars/pg+ws/2006/pem/ 

[25] Lukash V.N., Zh. Eks. Teor. iz 67,1594 (1974) 

[26] R. Conte, Physics Letters A 372, 2269 (2008) 



27 



